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Abstract 

Effects of the vacuum polarization leading to change of the effective grav- 
£^ | itational constant and particle creation in the early Universe are discussed. 

Gauss-Bonnet type coupling to the curvature is considered. Renormaliza- 
tion methods are generalized for such coupling and the iV-dimensional space- 
time. Calculations of creation of entropy and visible matter of the Universe 
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■^j- I by the gravity of dark matter are made. 
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O 1 Introduction 

o 

^ Quantum field theory in curved space-time is the generalization of well developed 

theory in Minkowski space-time, so it has some features of the standard theory 
together with new ones due to the curvature. It was actively developed in the 70- 
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ties of the last century (see books [UI2]) however some new results were obtained 
just recently and in this paper we'll concentrate on these new results as well as on 
some physical interpretation of the old ones. 

One of the main results of the theory for the early Friedmann Universe was the 
calculation of finite stress-energy tensor for particle creation. Let us begin from 
some remarks on what is meant by particle creation in curved space-time? In spite 
of the absence of the standard definition of the particle as the representation of 
the Poincare group in curved space-time where Poincare group is not a group of 
motions, one still has physically clear idea of the particle. The particle is under- 
stood as the classical point like object moving along the geodesic of the curved 
space-time. This classical particle however is the quasiclassical approximation of 
some quantum object, so the main mathematical problem is to find the answer to 
the question: what is the Fock quantization of the field giving just this answer for 
particles in this approximation? The answer to this question was given by us in 
the 70-ties for conformal massive scalar particles and spinor particles in Friedmann 
space-time by use of the metrical Hamiltonian diagonalization method. Creation 
of particles in the early Universe from vacuum by the gravitational field means 
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that for some small time close to singularity the stress-energy tensor has the ge- 
ometrical form, expressed through some combinations of the Riemann tensor and 
its derivatives while for the time larger than the Compton one it has the form of 
the dust of pointlike particles with mass and spin defined by the corresponding 
Poincare group representation. 

So our use of quantum field theory in curved space-time with its methods of 
regularization (dimensional regularization and Zeldovich-Starobinsky regulariza- 
tion) was totally justified by the obtained results. However still unclear was the 
situation with the minimally coupled scalar field where our method led to infinite 
results for the density of created particles, it is only recently that Yu. V.Pavlov [3 
could find the transformed Hamiltonian, diagonalization of which leads to finite 
results for this case. 

The effect of particle creation in the early Universe can have more general 
meaning concerning the origination of space-time itself [3]. It is reasonable to 
consider macroscopic Universe as classical approximation of the quantum Uni- 
verse. But what is quantum Universe? One can speculate that not only the metric 
but even the very existence of space-time points described by some coordinates 
is due to existence of operators of coordinates for some massive particles. In the 
"empty" quantum Universe there are no points! This is in some sense return to 
the Leibniz idea that space and time are relations between particles and don't 
exist without them. It is important to note [4J that space with space-like intervals 
between particles and time are necessary in order to measure or have information 
for Boolean - minded classical observer about observables described by noncom- 
muting operators. Noncommuting observables become commuting ones if one has 
some "copying" mechanism for the quantum system, when the same (identical) 
system is observed in different points of space separated by the space-like interval 
or at different moments of time. Time is needed for measuring noncommuting 
observables. This "copying" process makes possible contemplation of the quan- 
tum Universe with its complementary characteristics by the classical observer. 
"Copying" is the same as particle creation. That is why particle creation can be 
considered as origination of the classical Universe with its space and time from 
the quantum one. There is no necessity for quantization of gravity from this point 
of view, because space, time and metric are just some forms used by the classical 
observer dealing with the quantum Universe. Entropy of the Universe arising due 
to such "measurement" process is the other side of this process of particle creation 
and origination of space-time jlj. 

We use the system of units where h = c = 1. The signs of the curvature tensor 
and Ricci tensor are chosen such that = di Y l - k — T* 7 + r^r™ fc — T^T^, 
Rik = R\iki R — R\i where Y l - k are Christoffel symbols. 
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2 Some remarks on the physical interpretation 
of terms in vacuum polarization dependent on 
mass 

The calculations of the vacuum expectation value of the stress-energy tensor of 
the quantized conformal massive scalar, spinor and vector fields P3 |2] led to the 
expressions different for strong and weak external gravitational field. By strong 
field one understands the field with the curvature much larger than the one defined 
by the mass of the particle. For strong gravitational field it is basically polarization 
of vacuum terms which become zero if the gravitation is zero, for weak gravitation 
it is defined by particles created previously, so that if gravitation becomes zero 
particles still exist. Vacuum polarization for strong gravitation consists of three 
terms: the first is due to the conformal anomaly and it does not depend on mass 
of the particle at all, the second is due to the Casimir effect and is present for the 
closed Friedmann space, the third is vacuum polarization dependent on mass of 
the particle which is for scalar conformal particles |2] 

Here = Rn~ — Rgik/2 is the Einstein tensor, is the metrical tensor, 3? is some 
geometrical term of the dimension of the fourth degree of mass. It is interesting 
that the first term is the illustration of the Sakharov's idea [5] of gravitation 
as the vacuum polarization. If one takes the Planckean mass one has just the 
standard expression for the term in the Einstein equation. One can also think 
that gravitation is the manifestation of the vacuum polarization of all existing 
fields with different masses. However for weak gravitation for the time of evolution 
of the Friedmann Universe larger than the Compton one defined by the mass there 
is no Sakharov term, it is compensated. It is just manifestation of the fact, known 
also in quantum electrodynamics that vacuum polarization is different in strong 
and weak external fields being asymptotics of some general expression. So the 
physical meaning of this term in strong field is finite change of the gravitational 
constant, so that the new effective gravitational constant is different from that in 
the weak field, being 

I I m 2 
= 1 • 2) 

87rG e// 8ttG I44tt 2 1 ; 

So the larger is m, the smaller is G e ff. If m is macroscopic jH], then G e ff — > 
which means some kind of "asymptotic freedom" for gravity if one deals with 
macroscopic masses (like the mass of the Universe) quantum mechanically. 

The second term on the right side (0) describes what is now called "quint- 
essence" - cosmological constant, dependent on time. This "quintessence" is dif- 
ferent for different stages of the evolution of the Friedmann Universe. We calcu- 
lated it previously for the radiation dominated Universe, so that together with the 
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Sakharov term one has 



(T °) 



m 



48tt 2 U 2 



m 2 / a' 




167T 2 



4 

m 



ln(ma) + C + - + 




(3) 



Here a is the scale factor of the Friedmann space-time, r\ is the conformal time, 
C = 0.577 ... is the Euler constant. 

3 Creating of scalar particles in curved space- 
time 

One of unsolved problems of quantum theory in curved space is related to the 
definition of an elementary particle and vacuum in a curved space-time. The 
reason is that in the case of a curved space there is no group of symmetries similar 
to the Poincare group in the Minkowski space. If we believe that a particle is 
associated with a quantum of energy, then according to quantum mechanics, the 
observation of particles at an instant of time implies finding an eigenstate of the 
Hamiltonian. The Hamiltonian diagonalization method [2] takes this into account 
automatically. 

The use of the Hamiltonian constructed from the metrical energy-momentum 
tensor, first proposed by A. A. Grib and S.G.Mamayev [7] was successful in the 
homogeneous isotropic space-time for the conformal scalar fields. But such Hamil- 
tonian leads to the difficulties related to an infinite density of created quasiparticles 
in the nonconformal case The energy of such quasiparticles differs from the 
oscillator frequency of the wave equation jUJ QUI- The nonconformal case is im- 
portant for investigation by many reasons. Vector mesons and gravitons can be 
such particles. The additional nonconformal terms can be dominant in the vacuum 
expectation values of the energy-momentum tensor [TI]. In Ref. [2j the modified 
Hamiltonian was found so that the density of the particles corresponding to its 
diagonal form and created in the nonstationary homogeneous isotropic space-time 
is finite. In this section one gives the generalization of the method, proposed in 
Ref. for the cases of conformally static spaces and coupling with curvature of 
the general form (see Ref. 

We consider a complex scalar field tp(x) of the mass m with Lagrangian 




y[\g\ [g ik diip*d k ip - (m 2 + V g ) if* if 



(4) 



and corresponding equation of motion 




(5) 



where V« are the covariant derivatives in iV-dimensional space-time, gr = det(^), 
V g is the function of the invariant combinations and curvature tensor R 1 ^- 
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The equation (J3J is conformally invariant if m = and V g = £ C R, where £ c = 
(N — 2)/ [4 (iV — 1)] (conformal coupling). The case V g = is the minimal coupling. 
In this section we investigate the case of conformally static metric of the form 



ds 2 = dt 2 — a 2 {t) 7 Qj a(x) dx a dx 13 = a 2 {rf) (drj 2 — 7 aj g(x) dx a dx^ 



(6) 



where a, [3 — 1, . . . , N — 1. It is realized, in particular, for the homogeneous 
isotropic space-time. The equation (jSJ) for the functions coordi- 
nates (rj, x) takes the form 



(m 2 + V g 



N -2 



(2c + (N — 2)c 2 ) 



<p = 0, 



(7) 



where the prime denotes the derivative with respect to "conformal" time rj, c = 
a'/a , Atv_i = 7~ 1 / 2 d a ( v /7 7 a/3 <9 / 3), 7 = det(7 aj g). Let us assume that Laplace- 
Beltrami operator Ajv-i has the complete orthonormal set of the eigenfunctions 
$j(x) 

A A r_ 1 $ J (x) = -A 2 (J)$ J (x), (8) 

where J is set of iV — 1 indices (quantum numbers). The full set of the eq. (J7J) 
solutions can be found in the form <p(x) = gj(rj)<&j(x) , where 



g>>( V ) + n 2 (r ] )g J (r ] ) = 0. 



Q is the oscillator frequency: 



tf( v ) = (m 2 + V g ) a 2 - ^— - (2c' + (N- 2)c 2 ) + A 2 (J) . 
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(9) 



(10) 



The expression gj{rj) g*j{rj) — gj(v) 9j(v) is the first integral of eq. (JSJ). Further 
one normalizes the solutions of eq. © by condition 



9 J 9 j -9j9j = -2* 
For quantization, we expand the field <p(x) 



^(x) = / J) 



where dfi(J) is the measure on the space of Ajv_i eigenvalues, 



(12) 



= ^=gM $](x) , <fc>{x) = [0?>(x) 
and impose the usual commutation relations on the operators a/j , ctj . 



;(-), 



*(+)/ 



a (+) 



>JJ' 



a j , a j? 



a (±) a (±) 







(13) 



(14) 



We construct Hamiltonian as canonical one for variables <f(x) and ip*(x). If we 
add A^-divergence (dJ l /dx l ) to Lagrangian density (0J), where in the coordinate 
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system (77, x) the iV-vector (J 1 ) = (^j c(p* (p(N — 2)/2, 0, ... ,0), the motion 
equation (J7J) is invariant under this addition. By using the Lagrangian density 
L A (x) = L(x) + (dJ l /dx 1 ), we obtain by integration on hypersurface S : rj = const 



of Hamiltonian density h(x) 
modified Hamiltonian 



<p'(dL A )/(d<p') + <p*'(dL A )/(d<p*') - L A (x), the 



H(n) 



+ 



h(x) d N ~ 1 x 



(m 2 + V g ) a 2 



N -2 



(2c' + (JV — 2)c 2 ) 



99 



(15) 



Hamiltonian ()15|) is expressed in terms of the operators 



(±) *(±) 



by 



= fdf,(J)\E J (r ] )(a^a ( j- ) + a ( fV +) 



V J 



where 



EM 



\g'j\ 2 + ^\gj\ 



2 



9 j + ^ 9 j 



2„2 



(16) 
(17) 



here we used the normalization condition (jllj) and the special choice of the eigen- 
functions so that for an arbitrary J such J exists, that $}(x) = -$j$j(x). In 
accordance with condition of the completeness such set can be selected by the 
redefinition of eigenfunctions $j(x). The Hamiltonian diagonalization for initial 
instant t]q and normalization condition (|llj) give (under Q 2 (rj ) > 0): 



g'Avo) = i^(vo)gj(vo) , \gj(m)\ = ^ 1/2 M ■ 



(18) 



The Hamiltonian diagonalization for an arbitrary instant r\ is realized in terms of 

the operators b^, £>j^ related to the operators a^, via the time-dependent 
Bogoliubov transformations 



4~ } = ^M^i^-^mjbfivl V = aMb { iXv)-PM$jbf{r)l (19) 

where the functions aj, (3j satisfy the initial conditions |o; j(?7o) | = 1, /3j(vo) = 
and the identity |aj(r/)| 2 — \Pj{rj)\ 2 = 1. Substituting expansion (fTUj) in (fTo^) and 

requiring that the coefficients of the nondiagonal terms bj bj vanish, we obtain 



- 2ajf3j$jEj + a 2 Fj + (3 2 d 2 F* 



m 2 = ^(\9'.j\ 2 + tf\gj\ 2 )-\, (20) 



H{n) = J du.(J) nfo) 6^6^) • (21) 

Therefore, the energies of quasiparticles corresponding to the diagonal form of the 
Hamiltonian (JTfij) are equal to the oscillator frequency 

The vacuum state for the instant n is defined by bj '(t)) \0 n ) —bj 0. 
The state |0) = \0 VO ) contains \/3j(r})\ 2 quasiparticle pairs corresponding to the 
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operators b j (77) in every mode EB] • The density of created particles is propor- 
tional to Jdfj,(J) |/?j| 2 . The function S(rj) = \f3j{rj)\ 2 obeys the integral equation 

1 rv rvi 
S{v) = t; / dTuwfa) / dr] 2 w{n 2 ){l + 2S{r] 2 )) cos[2 6(^,^1)], (22) 

where ^ 

Using the first iteration and taking into account that 0(r/2,?7i) — * A (771 — 772) in 
A — > 00 one can represent (}2*2*|) as 



v 

w(t)i) exp(2zAr/!) (£771 

vo 



(23) 



Consequently, one can see, that S ~ A -6 in A (J) — ► 00. 

Therefore, in 4-dimensional space-time with metric the density of created 
particles is finite. So it is shown differently from [HI Q3], that even for minimal 
coupling the result is finite! From eq. (jUJ) for S(rj) one can see that under V g = (,R 
the densities of scalar particles created in Friedmann radiative-dominant Universe 
don't depend on £. So the results will be for t ^> m" 1 the same as for conformal 
coupled particles obtained earlier [UI2] ( see ec L- m this paper). 



4 Scalar field with Gauss-Bonnet type coupling 
to the curvature 

Usually for scalar field in curved space-time one writes the eq. (0) with V g = 
The condition £ = £ c consistent with conformal invariance of the equations for 
massless field is considered sometimes as preferable (see for example O ITB]). 
For scalar fields in inflation models JH] minimal coupling £ = usually is used. 
In general case one can take arbitrary values of £ [TT) ITT] . Renormalization of 
interacting fields in curved space-time is inconsistent with conformal invariance 
not only of the effective action but also of usual action ^3]- Models with arbitrary 
£ without conformal invariance on the classical level can be generalized by adding 
of the quadratic in curvature terms 

V g = + (R 2 + KRijW + xR^kiR ljkl + ■■■ (24) 

Nonzero coupling constants in these terms in (|24|) having the dimension (mass) -2 
lead usually to derivatives of the third and fourth order in the metrical energy- 
momentum tensor (EMT) and in Einstein equations. As it is known 18J such 
terms with higher derivatives even for small coefficients lead to the radical change 
of the theory. If one has the condition that metrical EMT of the scalar field does 
not contain derivatives of metric higher than the second order, one can take for V g 

V g = £R+(R* B , where Rq B = f Ri mpq R lmpg — 4Ri m R lm + R 2 . (25) 
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On even-dimensional spaces due to the Gauss-Bonnet theorem Euler charac- 
teristic of compact orientable manifold M with Riemann metric is equal to 



X(M) 



M 



E(x)^gd N x 



(26) 



where E(x) = -(4tt)- 1 j R(x) for N = 2 and E(x) = R^ B /(128n 2 ) for N = A. So 
the case with V g defined by (|23j) can be called Gauss-Bonnet type coupling. 

The condition of absence of higher derivatives of metric was considered ear- 
lier [20] as basic for multidimensional generalizations of the gravity theory. It was 
shown that the demands of symmetry in indices, covariant conservation and the 
condition that the tensor generalizing Einstein tensor Gi k is constructed from 
and only its first and second derivatives give the Einstein equations with cosmolog- 
ical constant in dimension four. But in dimension larger than four the correspond- 
ing equations can be obtained from the action which is a sum of contributions 
associated with continuation to this dimension of Euler characteristics of all lower 
even dimensions [20j. The demand of absence of ghosts in low energy approxima- 
tions of string theories gives the Einstein-Gauss-Bonnet gravity with Rq B density 
in action for higher dimensions [21]. The Gauss-Bonnet coupling of scalar field 
with gravity was also considered in dilaton theories (see for example [2*2"]). 

Taking variational derivatives of the action for scalar field with coupling 
one obtains the following expression for metrical EMT: 



where 



T, 



ik 



E. 



ik 



di(p*d k (p + d k <f*di<p - g ik d l (p*di(p + g ik m 2 ip*(p - 

2£ (G tk + V, V fc - g ik V l Vi) (<p*<p) - 2C (E ik + P ik ) , (27) 



5 R 



\g \ d N x 



GB 



p p Imp _ 9ik_ p ftlmpq _ 
-•-^ilmp-^k 2 ^Impq^ 



g\ 5g ik 

nk 

lmp ^C lmpq C lm ™ -(N- A)^H lk 



— 4R lm R lmk - 4RuR l k + 2g lk R lm R lm + 2RR lk - 9 fR 2 

— or 1 c 

— ^^ilmp^ k 



(28) 



% k 



ik 



i?V,V fc + 2^^' + 2g lk R lm V l V m - g lk RV l V 



l\-7m 



+ 4 



4i2 {( iVfc)V' - 2R ilkm V% 
N-3 



-AC tlk mVV m + 
R lk VN l + g lk R lm V l V m - 



N-2 



- 2%V fe) V< + 



NR 



2(N- 



AC ilkm R lm , 2{N- 3) 



(29) 



N-2 



(N-2Y 



2R ih R k 



NRR ik 
N - 1 



9ik RlmR 



Irn 



(N+2)R 2 
4(iV- 1) 



(30) 
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n D 2 / \ 2Rg l [ i g k ] m 

^iklm — tiiklm + Jj— ~ ytlm[i 9k]l ~ til\i9k\m J + ^ _ fJTj^Z 2) ' ^ ' 

parentheses (square brackets) in indices mean symmetrization (antisymmetriza- 
tion). Expressions for EMT using conformal Weyl tensor Ciki m are convenient for 
calculation as in conformally flat (Cn-im — 0), for example homogeneous isotropic 
spaces, as in Ricci flat (Rik = 0) spaces. For iV = 2, 3 one has Rq B = and 
there are no new effects due to C, ^ 0. In 4-dimensional space-time = (see 
Ref. |2HJ), but Pik(ip*ip) 7^ for general metric and (f(x) ^ const. 

Vacuum averages of EMT for quantized fields are divergent. For the analysis of 
the geometrical structure of divergences of vacuum averages of EMT one can use 
dimensionally regularized effective action. For complex scalar field with eq. (0) in 
loop approximation the effective action can be written |23J [2H] as 



S, 



eff 



J L eff (x)J\g~\d N x, (32) 



where 

L eff {x) 



(4?r) -iV/ 2 g ^ m N -2 jT ^ _ ^ (33) 



j=0 



a (x) = 1 , ai {x)= l -R-V g = ( j^^ +^jR- <R 2 gb , (34) 

a M _ Ri mpq R lmpq Ri m R lm R 2 V%R RV g V g 2 

a2{x) ~ 180 180 - + 72 ~ "10 6~ T _ Q~ ' (35) 

N is space-time dimension which is taken as some variable, analytically continu- 
ated in complex plane, e is some complex parameter, M is some constant of the 
dimension of mass necessary for correct dimension of L e tf, (length) ~ N ° for 
N = N — 2e, T(z) is gamma-function, A£ = £ c — £. Using (|3ip for Gauss-Bonnet 
coupling (J23j) one can write a 2 as 

/ (N-A)(N-Q) A^N-A) (AQ 2 \ 2 (jV-2)q mp(? C^ 
21 J ^ 480(^-1)2 12(iV-l) + 2 y + 240 (N - 3) 



720(iV-3) V 12 ^" 1 ) 7 2 6 V5 

First [iVo/2] + 1 terms in are excluded to get the renormalized L e ff ([b] 
denotes the integer part of the number b). By variation in of terms j = 0, 1, 2 
in effective action one obtains terms subtracted from vacuum EMT 
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TikA 1 } 



m 



N "- 2 /4ttM 



2\£ 



r 1 - 



X 



N 



x 



m N »- 2 / 4ttM 2 \ £ 
2^0-1^0/2 1 m 2 J 



ifc 



G 



ik 



4(E,, 



ik 



(N-2) V 3 ^- 1 ) (N-4) 



Tik,e [2] ; 



m 



N - 



(47T) 



, a V m 2 J |360(7V-3) 



(iV-6)S ifc + 3(JV-2)Wi fc 



+ 



r(4-f) r(s-f) 9 / ^ 

60(7V-1) 2 ? 3(N-1) { V 2 

r(3-f) 

3(iV-l) 



C 



+ 



c 



where 



ik 



Wi 



ik 



sJc lmpq C lm ™J\g~\d N x 



'\g\5g ik 

NRR ik , NR 2 g ik N-2 



— rr _l 4(^V3) / /m im 



(38) 



(39) 



+ < 2 r( 2 ^) 

A£ ^ (i^ + V, V fc - V, V') i^ fl + (£ ife + P, fc ) I , 



5[R 2 ^g~\d N x , , \ / 1 . 

' =— = 2N l VkR-g l kV l V l R)+2R(R lk --Rg l k), (40) 

5g tk ^ / V 4 



+ 



+ 



ViV k R + 



V'V^-V'Vji^fc . (41) 



2(iV-l) 8(iV-l) 2(iV-l) 2(7V-1) 
In conformally flat case Ciki m = 0, E^j '(4 — iV) = ^Hik, Wik = and so 
m N »~ 2 (AixM 2 \ S \ >^f_ N\„ r(3-f) / 



2^,-1^^0/2 I m 2 



(42) 
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Tik,e [2] 



UtcM 2 \ £ 



(47t)m,/ 2 y m 2 J 

N\f( N-A 



-r(4-f) 

90 (N - 3) 



+ 



A£ )R + (R GB 



+ 



r(4-f) r(s-f) , / N 

— V 1±_ + A t _V 2J_ / A .x 2 r ( 2 _ _ 

60(iV-l) 2 ? 3(iV-l) 1 I 2 



(43) 



-c 



+ ViV fc - & fe V z V ( i2£ B + P lk R 



r(3-f) 

3(JV-1) 

T 

Supposing that vacuum averages of EMT, i.e. (T^), are sources of the gravity 
field 

G ik + Ag ik = -8nG + ( T lfc ) V (44) 



C ( ^R G b + ^ ) #gb - A£ ( ( i? lfe + Vi V fc - </« V, V' ) Rqb + P ik R 



where A, G are cosmological and gravitational constants, T^ k is EMT of the back- 
ground, and taking into account eqs. (|3Tj) - (|3T)|) one comes to the following conclu- 
sion. First three subtractions from the vacuum EMT in iV-dimensional space-time 
correspond to renormalization of the cosmological and gravitational constants and 
parameters in quadratic, cubic and fourth degrees in curvature terms in the bare 
gravitational Lagrangian 



R-2A e 
16nG F 



+ a e R 2 GB + (3 £ R 2 + l£ C lmpq C lmpq + 5 £ RRq B + 6 £ R* B 



(45) 

Subtraction Tjfc j£ [0] due to (|37)1 . leads to infinite renormalization of the cosmo- 
logical constant A £ . From (|3*Hj) it follows that subtraction T^fl] correspond to 
renormalization of the gravitational constant G £ (finite for A" — > 4, £ = £ c ) and 
infinite renormalization of a £ (for ( ^ 0). Subtraction T^p] due to (|3*Bj) and ()39|) 
corresponds to renormalization of the parameters a e , @ e , j E , 5 e , 9 e . For £ = £ c and 
A" — > 4, the parameters /3 e , 5 £ have finite renormalization. 

Note that introduction of the term 0! £ Rq B in is consistent with dimen- 

sional regularization (E ik = only for integer AT = 2, 3, 4). 

For N -> 4 the products E ik V(l-(N/2)) and E ik V(2-(N/2)) have finite limits 
for arbitrary metric of space-time because the dependence of Ei k on A" is fraction- 
rational in analytic continuation on dimension, Ei k = for A^ = 4 and gamma- 
function has pole of the first order in points 0,-1. So the corresponding terms 
in (JHBJ), (jSUjl are finite and there is no necessity for subtracting them to get finite 
terms in (|44|) leading to back reaction of the quantized field in metric. However the 
effective action is divergent without such subtractions and the anomalous trace of 
vacuum EMT is different from the standard even for £ = 0. Due to the appearance 
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of such terms for different regularization procedures [21 it is convenient to have 
them in counterterms in vacuum EMT. 

Values of renormalized parameters are fixed by experiment and it can be [2] 
that renormalized parameters for non-Einstein terms in the gravity Lagrangian are 
zero. 

Now consider calculating of the renormalized vacuum EMT for iV-dimensional 
quasi-Euclidean space-time with metric 

ds 2 = dt 2 - a 2 it) dx 2 = a 2 (r?) (drj 2 - dx 2 ) . (46) 

In this case one can choose 

$j(x) = (2n)^ N ~ 1 ^ 2 e~ iXaXa , (47) 



J = {Ai, . . . , Atv-i} , —00 < \ a < +00 , A 



N-l 



\J2 X li 9j(v) = 9x(v) ■ 
\ 0=1 



It is convenient to express vacuum averages of the EMT operator for Fock vac- 
uum |0) annihilated by operator aj~' (|T4*|) through linear combinations of 
functions g x , g* x : 



s _ Wx\ 2 + n 2 \g x \ 2 1 ^IffAl 2 - 1 *^ 
AVI 2' 2Q 

which due to satisfy differential equations 



V 



d(g* x 9x) 
2dn 



(48) 



S'=^-U, U' = ^-(1 + 2S) -2QV , V' = 2ttU. 

i i L i U 



(49) 



Taking initial conditions S(r] ) = U(rj ) = V(r]o) = following from (fTH|) the 
eqs. (j4*9*|) can be written as integral equations of the Volterra type ()22|) and 



U{ri)+iV{rj)= / w(r ) i)(l + 2S(r ]1 ))exp[2ie(r] 1 ,r ] )}dr) 1 



(50) 



no 



Substituting expansion ^ in eq. (J2ZJ) and using (ffljl . (J2HI) (JHOl) , (HHJ), one obtains 
the following (diverging) expression for vacuum EMT averages 



(0|T ifc |0) 



B 



a N ~ 2 jo 



N 1 " n- \N-2j\ 

T ik A a\ 



(51) 



where B 



N : 



2^-3 7r (JV-D/2r((JV-l)/2) 



T 00 = n[S + -)+(N-l)[A£-(c 2 



cV+[c' + {N-2)c 2 )-[S+-U + - 



1 1 



(52) 
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T, 



a/3 



>a(3 



A 2 S + ^ 



x z ) — 



+ 



A£(N-1) - C (N+l)c 2 - 2d 



cV-2 A£-Cc 2 mU- (53) 



A£(N-l)d - ((N-l)c 2 (3c'-2c 



C = (a- 2 2(N-2)(N-3). 

The vacuum expectation (|5T|) has [N/2] + 1 different types of divergences: 
~ X N , X N - 2 ,...,\n\ if N is even, and ~ X N , X N ~ 2 , . . . , A if JV is odd. 
For renormalization we are using n-wave procedure proposed by Zeldovich and 
Starobinsky [27\ 



(0|T ifc |0), 



B 



N 



,N-2 



N-2 



1 d l 
Tik[l] = 77 hm 



[N/2] 

Tik - ^2 Tik [l] 
1=0 

-Tik(nX,nm)) 
n / 



dX . 



(54) 



(55) 



To obtain the explicit expression r^fZ] write the series for S,U,V in reverse 
degrees of n for A — > nX, m — > nm, n — > oo: S = Y^k=i n ~ k Sk, ■ ■ ■ Taking 
subsequent iterations in integral eqs. (J22J, (|5Uj) and using the stationary phase 
method one obtains for first different from zero terms 



V 1 = W, U 2 = DW , S 2 = -W 2 

4 



-It - D-W- -d(^- 

uj 6 



04 



- W DW — D W — D \^—D ^-^ ] j 

3 



q 



2u 2 



:DW. 
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TU 4 + \ (DW) 2 - - WD 2 W - - ujWD (^) 
4 2 4 Vu; 3 / 



where 



g =(A^i?-Ci?Gs) a2 ; w = (mV + A 2 ) 1/2 , W 



2^ 



D = — — 

2u drj 



(56) 

(57) 
(58) 

(59) 



Note that in (j5T?j) - (p)8*j) terms nonlocal on time (dependent on n and rj ) are ex- 
cluded. These terms are absent if Vi(rj ) = V^o) = ^2(^0) = ^4(^0) — which is 
supposed further. In particular, nonlocal terms are absent if first 2[N/2] derivatives 
of the scale factor a(rj) of metric are zero in the initial moment of time. 
Using (|52 jl . (f53 1l . ©-(0 for r ik [l] one obtains 



roo[0] 



CO 

2 



T a p[0] 



X 2 



2{N-1)uj 



r 00 [l] = u;S 2 + (N-l) (AC - (c 2 ) cVi + 



N-l 



A£2(N—2)c 2 + ((4-3N)c 4 



(60) 



(61) 
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W 1 ] 



Ja/3 



(N-I)lu 



+ 



+ 



+ 





1 


+ 






] = 


< 


N- 


-l r 



Af(JV-l) - ((^N+l)c 2 - 2d J cVi - 2 (A£ - Cc 2 ) wtf 2 + (62) 
A£(iV-2) (c 2 - 2c') + Cc 2 (3iV-4) (2d - \c 2 



Alu 2 



U 2 + ^) + {N-l)(Ai-Cc 2 )cV,+ 



1 



N-l 



A£2{N-2)c 2 - C(3iV-4)c 4 



w P 4 + 4^ + 16c^ / 



m 2 a 2 



+ 



AC(iV-l)-C (iV+l)c 2 -2c' 



cv:. 



2w 

A^ - Cc 2 



4w 4 



2wf7 4 



A^(iV-2)(c 2 -2c') +C(3A^-4)c 2 (2c'- -c 2 



(63) 

+ 

+ (64) 



These expressions give all information on subtractions for N = 4, 5. New coun- 
terterms appear for N > Q. For conformal scalar field they are given in Ref. j2H|- 
The renormalized due to (p>4*j) vacuum EMT is covariantly conserved. This is seen 
from equations V l (T ik /a N ~ 2 ) = and V l (Tik[l)/a N ~ 2 ) = following from 

052), (El, (EOD (EH)- 

To see the geometrical structure of the n-wave procedure counterterms let us 
make as in [29J the dimensional regularization. For calculation of integrals in 
dimensionally regularized counterterms 

TikAQ = 4k( M ) 2e r* N ~* T *M dx > ( 6s ) 



a 1 " " jo 
with T ik>£ [l] defined by pijl - (jH4) with putting N 



00 r(t±l)Y(n 
x k (l + x i)-Pdx = - [ 2 } {P 



N — 2e, one uses the equality 

(66) 



fc+i 
2 



10 2T(p) 

In cases when the integral in the left hand side of (|66|) does not exist in usual sense 
it is taken as the analytic continuation of the right side of (|66|) for corresponding 
values of k and p. As the result of calculations one obtains expressions (pT7J). (J4*2~j) . 
(J4*3~j) for the 0-th, 1-st and 2-nd counterterms n-wave procedure. So the geometrical 
structure of first three subtractions in the n-wave procedure and in the effective 
action method occurs to be the same. 

In conclusion note that Gauss-Bonnet coupling terms Rqb can pl & y important 
role in the early Universe. Effects of the nonzero value of ( for scalar field can be 
observed in black hole radiation, in parameters of the so called bosonic stars [3*Uj . 
Explicit value of ( can be checked by experiment. 
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5 Superheavy particles in the early Universe 

It is known [21 EH1 E21 that the number of particles with mass of the order of the 
Grand Unification scale created by gravitation in the early Universe described by 
the radiation dominated Friedmann metric is of the Dirac-Eddington order, i.e. 
of the observable order for the visible mass. On the other side it is clear that if 
superheavy particles after their creation continued to be stable for large enough 
time they will lead to the collapse of the Universe governed by the radiation dom- 
inated metric in the short on the cosmological scale time for closed Friedmann 
space or lead to the unrealistic scale factor for the open space. So the idea was 
proposed that these superheavy particles must decay on quarks and leptons with 
CP- noninvariance leading to the observable baryon charge of the Universe before 
the time when the energy density of the created superheavy particles will become 
equal to that creating the background metric. If superheavy particles have nonzero 
baryon charge then their decay in analogy with decay of neutral i^-mesons will 
go as decay of some short living and long living components. Supposing that the 
lifetime of long living components is of the cosmological order but their number 
was diminished in comparison with the number of the short living components 
due to their interaction with the baryon charge created previously similar to the 
well known regeneration mechanism for i^-mesons one can speculate about their 
existence today as cold dark matter. Rare events of their decays can be identified 
as experimental observations of high energetic cosmic rays [33] with the energy 
higher than the Greizen-Zatsepin-Kuzmin limit j^l]. Here we shall discuss differ- 
ent possibilities of the role of superheavy particles with the mass of the Grand 
Unification scale in the early Universe. 

1) It is natural to think that some inflation era took place before the Friedmann 
stage. Some inflaton field probably manifesting itself as the quintessence in the 
modern epoch after the quasi de Sitter stage led to the dust like or to the radiation 
dominated Friedmann Universe. Usually it is supposed that the inflaton field does 
not interact with ordinary particles and can be some manifestation of the non 
Einsteinian gravity for example due to high order corrections. So even if it decayed 
on some light "inflaton" particles the primordial inflaton field can form hot dark 
matter but not the visible matter and entropy present in background radiation. 
Our idea is that inflaton field was the source of Friedmann metric with some 
small inhomogeneities, but visible matter and the entropy of the Universe were 
created not by the inflaton field itself but by the gravitation of this inflaton field. 
That gravitation created pairs of superheavy particles. Short living components 
decayed in time of the Grand Unification scale and led to the nonzero baryon charge 
observed today as visible matter. If long living components had the lifetime of the 
order of the "early recombination era" then the energy density of created long 
living particles soon became equal to that of the background inflaton field (hot 
dark matter). Then the decay of all long living components led to the observable 
entropy of the Universe. Here it is supposed that the energy density of the inflaton 
field led to the observed cosmological scale factor, so it is evident that the created 
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entropy due to our mechanism will be of the observable order. 

2) The other possibility is to put the hypothesis discussed by us earlier 
that not all long living components decayed and formed the entropy but some 
part of them survived up to modern time as cold dark matter and superheavy 
particles are observed in cosmic rays events. Then it is natural to suppose that 
the lifetime of the long living component is of the cosmological order but the large 
part of them regenerated into short living components due to interaction with the 
baryon charge in time shorter or equal to that of the "early recombination era" 
and entropy appeared due to this decay. 

Now let us give some numerical estimates. Total number of massive particles 
created in Friedmann radiation dominated Universe (scale factor a(t) = dot 1 / 2 ) 
inside the horizon is as it is known QjJ |2| : 

N = n^ s \t) a 3 {t) = b (s) M 3/2 a 3 , (67) 

where ~ 5.3- 10~ 4 for scalar and &( 4 / 2 ) ?s 3.9- 10~ 3 for spinor particles (N ~ 10 80 
for M ~ 10 14 GeV, see Ref. 0). Radiation dominance in the end of inflation era 
dark matter is important for our calculations. If it is dust like the results will be 
different (see further). For the time t 3> M~ l there is an era of going from the 
radiation dominated model to the dust model of superheavy particles 

If M ~ 10 14 GeV, tx ~ 10~ 15 s for scalar and tx ~ 10~ 17 s for spinor particles. 
Let us call tx - "early recombination era". For dust like ending of inflation era 
one has N ~ M (see Ref. p6 ) and therefore the ratio of the X-particles energy 
density ex to the critical density e cr u does not depend on time (ex < £ C rit for 
M < Mpi). 

Let us define d — the permitted part of long living X-particles — from the 
condition: on the moment of recombination t rec in the observable Universe one has 
de x (t rec ) = e CTit (t Tec ) . It leads to 

M Pl \ 2 1 



d 64vr^) { M ) s/MUTc (69) 

For M = 10 13 -10 14 GeVone has d « l(T 12 -l(r 14 for scalar and d w l(T 13 -l(r 15 
for spinor particles. So the life time of main part or all X-particles must be smaller 
or equal than tx- 

Now let us construct the model which can give: a) short living X-particles 
decay in time r q < tx (more wishful is T q ~ tc ~ 10~ 38 — 10~ 35 s, i.e. Compton 
time for X-particles) b) long living particles decay with 77 w tx- Baryon charge 
nonconservation with CP-nonconservation in full analogy with the i^°-meson the- 
ory with nonconserved hypercharge and CP-nonconservation leads to the effective 
Hamiltonian of the decaying X, X - particles with nonhermitean matrix. 
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For the matrix of the effective Hamiltonian H = {Hij}, i,j = 1,2 let H\\ =i?22 
due to CPT-invariance. Denote e = (V-^12 — y/Hn ) / (V-^12 + V H21 ) • The 
eigenvalues Ai,2 and eigenvectors | ^1,2) of matrix H are 



Ai, 



In particular 



H 



l*i,2) 



/ E 



1 



ifi2 + Hoi 1—6 



2(1 + kl 



2 1 + e 2 ' 
[(l + e)|l>± (l-e)|2)]. 



(70) 
(71) 



4 V T 9 1 



l+£ 
l-£ 



- r, 



\ l+e 



.-1 



_-l 



(72) 



/ 



Then the state \^\) describes short living particles X q with the life time r q 
and mass E + A. The state ^2) is the state of long living particles Xi with life 
time 77 and mass E — A. Here A is the arbitrary parameter — E < A < E and it 
can be zero, E = M. 

So for the scenario 1) it is sufficient to take 77 ~ tx- 

In scenario 2) the small d ~ 10~ 15 — 10 -12 part of long living X-particles with 
Ti > tjj m 4.3 • 10 17 s is the age of the Universe) is forming the dark matter. 
The decay of these superheavy particles in modern epoch can give observed ultra 
high energy cosmic rays. Using the estimate for the velocity of change of the 
concentration of long living superheavy particles [^7] \h x \ ~ 10~ 42 cm~ 3 s _1 , and 
taking the life time 77 of long living particles as 2 • 10 22 s, we obtain concentration 



-20 , 



at the modern epoch, corresponding to the critical density for 



rix ~ 2 • 10 cm 
M = 10 14 GeV . 

Let us use the model with effective Hamiltonian (|72j) where 77 > tjj and take 
into account transformations of the long living component into the short living 
one due to the presence of baryon substance created by decays of the short living 
particles in analogy with the regeneration mechanism for i^°-mesons. 

Let us investigate the model with the interaction which in the basis |1), |2) is 
described by the matrix 

/0 ° ^ (73) 







-Z7 



The eigenvalues of the Hamiltonian H + H d are 



Af, 2 = B-J(r- + rr' 




A--W. 



T 



(74) 



In case when 7 <C r 1 for the long living component one obtains 



E-A-lrr'-iZ 



(75) 
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ta-t ft 



7(t) dt 



(76) 



|*2(t)|| 2 = ||* 2 (to)|| 2 exp 

in Jt 

The parameter 7, describing the interaction with the substance of the baryon 
medium, is evidently dependent on its state and concentration of particles in it. For 
approximate evaluations take this parameter as proportional to the concentration 
of particles: 7 = an (0) (f). Putting n = 2 ■ 10 22 s, t < t v , a(t) = a Vi by (JHZJ) 
one obtains 



1*2(011 = ll*2(*o) II 2 exp 




(77) 



So the decay of the long living component due to this mechanism takes place close 
to the time to- One can think that this interaction of X/ with baryon charge is 
effective for times, when the baryon charge becomes strictly conserved, i.e. we take 
the time larger or equal to the electroweak time scale, defined by the temperature 
of the products of decay of X q . This temperature is defined from Mn^{r q ) « a~T 4 
and is given by 



T(t)= Uw) -kjT' (78) 

where ks is Boltzmann constant, Ni is defined by the number of boson Nb and 
fermion Np degrees of freedom of all kinds of light particles: iVj = Nb + |Xp (see 
Ref. [38). At time tx this temperature is equal to 

If r q = 1/M and Ni ~ 10 2 - 10 4 , then for spinor X-particles T(t x ) w 300 - 
100 GeV, i.e. the electroweak scale for created particles (which is however different 
from that for the background). 

So let us suppose to « tx- If d — is the part of long living particles surviving up 
to the time t {tjj > t ^> t c ) then from ([69)1 and (J77J) one obtains the evaluation 
for the parameter a 

-3 hid M 2 Pl 
a ~ 128tt(^)) 2 M 4 " ' ( } 

For M = 10 14 GeV and d = 10" 15 one obtains a w 10" 30 cm 3 /s. If r q ~ 10" 38 - 
10~ 35 s then the condition j(t) <^ t^ 1 used in Eq. (f?5|) is valid for t > tx- For this 
value a we have 7(%) ~ 10 -36 s _1 . So one can neglect this mechanism for 

the decay of the long living component of X-particles for the modern epoch while 
for early universe at to ~ tx it was important. The observable entropy in this 
scenario is created due to decay of Xi on quarks and antiquarks at the time tx 
when the Grand Unification symmetry is totally broken. Baryon charge is created 
at t q which can be equal to Compton time for X-particles tc ~ 10~ 38 — 10~ 35 s. 

Our scheme is the same for the scalar particles and the fermions. The super- 
heavy fermions are used, for example, in some models of neutrino mass generation 
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(the see-saw mechanism) in Grand Unification theories [351 EH]- New experiments 
on high energetic particles in cosmic rays surely will give us more information on 
their structure and origin. 
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